Complex Probabilities on R as Real Probabilities on C and an Application to Path 

Integrals 

Don WeingartenQ 
IBM Research, P.O. Box 218, Yorktown Heights, NY 10598 

We establish a necessary and sufficient condition for averages over complex valued weight functions 
on R N to be represented as statistical averages over real, non-negative probability weights on C . 
Using this result, we show that many path- integrals for time-ordered expectation values of bosonic 
degrees of freedom in real-valued time can be expressed as statistical averages over ensembles of 
paths with complex-valued coordinates, and then speculate on possible consequences of this result 
for the relation between quantum and classical mechanics. 



For a wide range of quantum mechanical systems, if 
real-valued time is analytically continued to imaginary 
values, the weight assigned to each path in the path inte- 
gral formula for time-ordered vacuum expectation values 
becomes real and non-negative, permitting interpretation 
as a probability density. Paths distributed according to 
this probability density can then be generated by the 
Metropolis method, one version of which makes use of 
the Langevin equation. 

For the Langevin equation applied to paths for imagi- 
nary time path-integrals, however, the coefficients enter- 
ing the equation are themselves analytic functions of the 
action occurring in the path weight. As a consequence 
of this analyticity, the Langevin equation can be ana- 
lytically continued, at least formally, to complex-valued 
weights jl]] , among which are those of the original path- 
integral with real- valued time. For complex path weights, 
the coefficients in the analytically continued Langevin 
equation and the coordinates in the path ensembles the 
equation generates are complex. 

Numerical tests and some analytic results concerning 
the Langevin equation with complex coefficients yield a 
complicated picture For some systems the method 
gives accurate answers, for others values of the coordi- 
nates blow up as the equation's running time is increased, 
and in still other applications averages over complex 
Langevin trajectories converge but to wrong answers. 
A question which follows from these results is whether 
the failure of the Langevin equation for many complex- 
valued path weights is simply a limitation of the Langevin 
equation or, alternatively, arises from the nonexistence of 
equivalent real, non-negative weights on complex-valued 
trajectories. More generally, with no reference to the 
Langevin equation, under what circumstances can aver- 
ages over complex weights on real coordinates be realized 
as statistical averages over real, non-negative weights on 
complex coordinates? The main result addressing this 
problem so far is a proof |J that an average over real co- 
ordinates of a complex polynomial multiplied by a com- 
plex Gaussian can be realized as an average over a real, 



non-negative weight on complex coordinates. 

In the present article, by an explicit construction un- 
related to the Langevin equation, we establish a neces- 
sary and sufficient condition for averages on R N over 
complex weight functions, with no explicit restriction on 
functional form, to be represented as averages over real, 
non-negative probability densities on C N . For the par- 
ticular case of N = 1, we show this condition is fulfilled 
for all complex weights normalized to a total integral of 
1. 

As an application of our result for general N, we 
show that real-time path integrals for many systems with 
bosonic degrees of freedom can be represented as averages 
over ordinary statistical ensembles of paths through com- 
plex space. Among the systems which can be represented 
in this way are scalar field theories, in any space-time di- 
mension, with self-interaction cf> p , p even, field theories, 
in any dimension, with bosonic degress of freedom which 
range only over a bounded region, and therefore lattice 
gauge field theories over any compact gauge group. We 
then briefly speculate on possible consequences of this re- 
sult for the relation between quantum and classical me- 
chanics. 

The systems we consider all have time reduced to a 
finite, discrete lattice of points. We do not examine the 
existence of limits as the time lattice spacing is taken 
to zero. We also do not discuss algorithms for actually 
generating averages over the probability densities on C N 
which we prove exist. 

For c(xi, . . . , Xn) a complex- valued weight function on 
(xi, . . . , xn) € R N , t{z\, . . . , zn) a real, non-negative 
weight on (z\, . . . , z^) € C , and p(r\, . . . , rjv) a real, 
non-negative weight on [r\, . . . , rjv) € [0, oo) , define the 
expectation values 
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f(zi,...,z N )t(zi,...,z N ), (lb) 



< / > P = 
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/(n, ■ ■ ■ ,r N )p(r 1 , . . . ,r N ) (lc) 

with zj given by Tj exp(i9j) in Eq. [ll]. We require c(x), 
i(z) and p(r) to have non-zero total integrals, then choose 
Z c , Z t and Z p to give < 1 > c , < 1 >t, < 1 > p all the 
value 1. 

We will show that for any c(xi, . . . ,xn), a neces- 
sary and sufficient condition for the existence of a 
t(z\, . . . , zjv) with the same moments 
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for all integer rrij > 0, is that there exist a p(ri, 
satisfying the bound 
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mi 
1 



< <r™ 1 



'JV 



(2) 
,tn) 

(3) 



for all integer > 0. 

For notational simplicity, we consider first a detailed 
proof for JV = 1, then briefly summarize the generaliza- 
tion to arbitrary JV. 

To show that Eq. ^ implies Eq. ^ for some p(r), we 
take absolute values of both sides of Eq. p| and obtain 

\Z~ l / x m c(x)dx\ < Z^ 1 I r m p{r)dr, (4a) 
> Jo 

p(r) = r det(z). (4b) 



But Z p is identical to Z t . Thus Eqs. [| implies Eq. ||. 

We now assume Eq. |^ and construct a t(z) which sat- 
isfies Eq. |^. For z given by r exp(i9) as before and a 
real parameter A > 1, Eq. implies we can define the 
functions 



v-^ exp(-im8) < x m > c 

q\{V) = > — = !• 



m>0 



X m < r m >, 



s x (6) = {l + 2Re[q x {9)}}. 
Eqs. |B| imply for m > 
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— exp(im9) s x (9) = 
2n 
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With the definitions 



< / >xt = Z£ 



2ts 



d9 / drrf(z)t x (z), (7a) 
o Jo 

t x (z) = r- 1 p(X- 1 r)s x (9), (7b) 



and Z\ t chosen to make < 1 > xt — 1, Eqs. || and then 
yield 



< z m >xt 



< x" 1 > c 
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which is Eq. g for N = 1 . 

By a choice of A we can insure that t x (z) is non- 
negative. From Eq. [n] it follows that t x (z) is non- 
negative if s x (8) is non-negative. Eqs. || and || then 
yield 



s x (9) > l-2\q x (e)\, 

I < x n > c 



> 1 - 2 



\ m < r n > v ' 

m>l ^ 
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Eq. |9c| imples that for A > 3, t\{z) is non-negative. 

Consider now arbitrary N . The proof that Eq. ^ im- 
plies Eq. H for some p(ri, . . . , rjv) for arbitrary JV is an 
immediate extension of the JV = 1 proof. The proof that 
Eq. U implies Eq. ^ requires very little more effort. For 
general JV, q x {9) of Eq. pa] becomes 



q x (9\, . . . , 9n) — 

exp(-i J2i rriiOi) < x™ 



E 



'-JV 



< r l • • ■ r JV >P 
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Eq. |9c|, insuring the positivity of s x (z\, . . . , zjv) and 
therefore of Ja(-Zi, ■ ■ ■ , zjv) becomes 



ax{6t,...,e N ) > 3- 



2A 



JV 



(11) 
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(A-l) w ' 
the right side of which is positive for 

Returning to the case of JV = 1, we now show that a 
p(r) exists fulfilling Eq. ||for any c(x) with < 1 > c of lj|]. 
To do this, we make use of a theorem giving conditions 
for the existence of a p(r) such that for all non-negative 
integers m 

< r m > p = s m , (13) 

for a specified sequence of real numbers sq, s\, . . .. From 
the s m , for any non-negative integer n, define the (n + 
1) x (n + 1) matrices 



H 



2/i 



Tr2n+1 _ 

n ij — s i+i+l) 



(14a) 
(14b) 



with < i, j < n. Then according to a result of Stieltjes 
H a p(r) obeying Eq. |l^ exists if for all non-negative n 



det(H n ) > 0. 



(15) 
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For any c(x), we will construct a sequence so, s i,---, 
obeying both Eqs. || and [l^. Choose s to be 1. Now 
suppose so, s%, . . . , Sfe, have been found obeying Eqs. || 
and [[5] for all n < k. By isolating the dependence of 



det(H k+1 ) on 



c fc+i 



it is not hard to show that a pos- 



itive real bound b exists, such that Eq. |15| is fulfullcd 
for n = k + 1 for any s^+i > b. We choose Sk+i to be 
the greater of | < x k+1 > c | and b. By induction, the 
resulting sequence So, s±, . . . obeys Eqs. || a ndll5| for all 
m. It follows that a p(r) exists fulfilling Eq. |l3| and thus 
also Eq. ||. 

The result we have just proved can not be extended 
automatically to N > 1 since, for N > 1, the analog of 
Eq. [l5| is not a sufficient condition for the existence of a 
p(r%, . . . , tn) fulfilling the analog of Eq. [[3] fl. 

For each c(xi, . . . , Xn) obeying Eq. pf for some 
p(r%, . . . , rjv), the non-negative, real t\(zx, . . . , zn) we 
construct satisfying Eq. ^ is only one representative of a 
large family of such functions fulfilling Eq. El At the very 
least, given any p(r\, . . . , r^r) obeying Eq. an infinite 
family of distinct p'(r 1; . . . , r N ) which also obey Eq. || 
can be constructed by shifting some part of the weight 
carried by p(ri, ■ ■ . , rjy) out to larger values of 7*1 , . . . , rjy. 
For example, for positiveai, . . . ,ajv, define p'(r%, . . . ,rjv) 
to be zero except on the region r\ > ai, . . . ,Tn > o-n, 
where it is given by p{r\ — a\, ...,rjf — apj). Each mo- 
ment of p'(r\, . . . , rjsr) is greater than the corresponding 
moment of p(ri, . . . ,rjv), and thus satisfies Eq. [| By 
means of Eqs. |l^ - [l^, p'(ri, . . . ,rjy) can therefore be 
turned into a new t'C[zi, ■ ■ ■ ,zm) with moments satisfy- 
ing Eq. §. 

An expression for q\{6x, . . . , 6jy) somewhat simpler 
than Eq. [To] can be found if p(r±, . . . in Eq. ^ is 
a product pyrx) . . .p(rjy) with moments < r 
for any R > and integer m > 0, obey 



>p which, 



<r m >p > b R R r 



(16) 



for some bn > independent of m. Then with the defi- 
nition 



fx(y) = E 



y" 



m>0 H 



(17) 



Eq. takes the form 



< fx[xiexp(~idi)} . . . fx[x N exp(-i9 N )} > c -1. (18) 

For many quantum mechanical systems, a finite vol- 
ume, lattice approximation to the path integral in real- 
valued time leads to an expectation value fulfilling Eq. 
^ and therefore respresentable as an ordinary statisti- 
cal average over an ensemble of random paths through 
complex space. Among the systems which obey Eq. || 
for some p(x±, . . . ,xjv) are scalar field theories, in any 
space-time dimension, with self-interaction (f> p , p even, 
field theories, in any dimension, with bosonic degrees of 
freedom which range only over a bounded region, and 



therefore lattice gauge field theories over any compact 
gauge group. In addition, of course, for any system for 
which the real-valued time path integral is correctly han- 
dled by the complex Langevin equation, the probability 
distribution t(z\, . . . , zjv) given by the limit as running 
time becomes large of the ensemble of points in C N struck 
by the Langevin trajectory fulfills Eq. |[ Since Eq. |^ 



such a system then has a p(r±, . . . , r^) 



implies Eq. 
fulfilling Eq. 

As an example of a proof of Eq. |3j, consider the path 
integral for a single anharmonic oscillator with time re- 
duced to a periodic lattice of points labeled by a posi- 
tive integer t < T. With xt the position at time t, the 
time-ordered expectation of a polynomial f(xi, . . . ,xt) 
becomes < / > c of Eq. [la] with 




KX f 6 



(19) 



where fi is the particle's mass, k the anharmonic spring 
constant, and <5 the lattice spacing. Rotating each inte- 
gral over real Xt to an integral over exp(—iTr/8)yt with y t 
real, taking the absolute value of both sides of the rotated 
version of Eq. [TaL and using the inequality 



(yt-yt+1) 2 < 2(y? + yf +1 ), 



(20) 



we obtain 



N 



>d ^ w-M <r ™ t> P' ( 21a ) 



p(r) = exp 



V2r 2 
fid 



Kr 4 5 ,(21b) 



with < . . . >p given by Eq. |k] for N = 1. Eq. 21a then 
gives Eq. [}] with 



P{xi, 



,x T ) = 



A = Max[l, Zp 



(22a) 
(22b) 



The proof of Eq. g for an anharmonic oscillator ex- 
tends easily to a <f) p field theory, p even, for any dimen- 
sion of space-time and to bosonic field theories, in any di- 
mension, with with degress of freedom bounded by some 
positive real B. In a pr oof o f Eq. ^ for bounded degrees 
of freedom, p(r) in Eq. 21b becomes 5(r — B). 

Finally, we speculate briefly on the possibility that the 
degrees of freedom of the real world might actually be 
complex- valued, distributed according to a probability 
which, as we have shown possible in many cases, yields 
averages agreeing with those of conventional quantum 
mechanics. A version of this idea suggested by the com- 
plex Langevin equation is considered in Ref . |J7| . A prob- 
lem encountered immediately by this proposal is that 



4 



the macroscopic world, governed nearly by classical me- 
chanics, exhibits only real-valued positions. This prob- 
lem would be solved if, within the statistical ensemble of 
paths through complex-space, the paths for macroscopic 
variables happened nearly to obey classical mechanics 
and to have only small imaginary parts. A formulation 
of quantum mechanics incorporating this feature would, 
potentially, be capable of resolving many of the puzzles 
in the interpretation of quantum mechanics ^ ^ . A dis- 
cussion of the proposed resolution of these problems by 
decoherence criteria llO|, fXTl and of the difficulties which 



where Sk is 



these encounter appears in Rcf. 12 



The speculation we offer here, in effect, is that quan- 
tum mechanics is a version of ordinary statistical mechan- 
ics but for paths through a complex-valued coordinate 
space. Although this prospect may sound implausible, 
the harmonic oscillator provides a simple model of how 
such behavior might conceivably come about. 

Assume periodic boundary conditions for time period 
T = 2M + 1 and Fourier transform the oscillator coor- 
dinate x t at integer time t according to the convention 



x t 



E 

0<fc<M 



2irkt 2nkt 
a k cos{^^) + b k sin{—--) 



(23) 



with real coefficients a k and b k , b identically 
0. The time-ordered expectation of a polynomial 
/(a , . . . , a M , h,..., b M ) is < / > c of Eq. [la| for 



c(a , . . . , a M , h, . . . , b M ) 



exp 



E (rt+rt, 



Sk 



0<k<M 



(24) 



Sk = - cos(k)] - k5, 



(25) 



and S is the time lattice spacing. 

In Eq. [la| expressed as integrals over real a k and bk, 
for Sk positive the integrals can be rotated, respectively, 
to \flc k and \fidk with real and d k . For s k negative, 
the Ofc and bk integrals can be rotated, respectively, to 
^f—iCk and sf^ldk- If [i, however, is given a microscopic 
positive imaginary part e and a correctly chosen real part, 
there will be a single k for which s k has no real part and 
a microscopic positive imaginary part. For this k the 
integrals over and bk can be left pure real. Inverting 
the Fourier transform of Eq. a statistical ensemble 
of complex-valued trajectories xt results with real, non- 
negative probability weight. A typical trajectory in this 
ensemble consists of the sum a real part, with amplitude 
of order e -1 / 2 , oscillating at the frequency predicted by 
classical mechanics, and a complex part independent of 
e. 

The preceding construction is, effectively, an appli- 
cation to the harmonic oscillator path integral of the 
method of steepest descent. A similar construction can 
be carried out for a free scalar field theory in any di- 
mension of space time. For more complicated field theo- 
ries, the method of steepest descent can also be applied 
to obtain a t(x\, . . . , xn) which includes real-valued tra- 
jectories obeying the classical equations of motion, but 
the t(x\, . . . ,xn) produced in this way will only approx- 
imately fulfill Eq. g. Whether for interacting systems 
t{x\, . . . , xn) can be found which fulfill Eq. and in- 
clude real-valued trajectories obeying classical equations 
of motion is an open question. 
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